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Abstract: Precision tests of the Standard Model (SM) require not only accurate experi¬ 
ments, but also precise and reliable theoretical predictions. Triple vector boson production 
provides a unique opportunity to investigate the quartic gauge couplings and check the va¬ 
lidity of the gauge principle in the SM. Since the tree-level predictions alone are inadequate 
to meet this demand, the next-to-leading order (NLO) calculation becomes compulsory. In 
this paper, we calculate the NLO QCD + NLO electroweak (EW) corrections to the W±ZZ 
productions with subsequent leptonic decays at the 14 TeV LHC by adopting an improved 
narrow width approximation which takes into account the off-shell contributions and spin 
correlations from the W - and Z-boson leptonic decays. The NLO QCD+EW corrected 
integrated cross sections for the W±ZZ productions and some kinematic distributions of 
final products are provided. The results show that both the NLO QCD and NLO EW 
corrections are significant. In the jet-veto event selection scheme with Pj^j et = 50 GeV, the 
NLO QCD+EW relative corrections to the integrated cross section are 20.5% and 31.1%, 
while the genuine NLO EW relative corrections are —5.42% and —4.58%, for the W + ZZ 
and W~ ZZ productions, respectively. We also investigate the theoretical dependence of 
the integrated cross section on the factorization/renormalization scale, and find that the 
scale uncertainty is underestimated at the LO due to the fact that the strong coupling a s 
is not involved in the LO matrix elements. 
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1. Introduction 

After the discovery of the 126 GeV Higgs boson by the ATLAS and CMS collaborations 
[1,2], the main tasks of further experiments at the Large Hadron Collider (LHC) are to 
determine the Higgs properties, test the predictions of the Standard Model (SM), and 
search for new physics beyond the SM. One of the important experiments for testing the 
SM is to measure the gauge couplings in the SM and check the validity of the gauge 
principle. The multiple gauge boson productions at the LHC can be used to determine 
the gauge boson self-couplings and help us for better understanding the electroweak (EW) 
symmetry breaking. The theoretical predictions for most multiple gauge boson productions 
at the LHC have been computed up to the QCD next-to-leading order (NLO) so far. The 
NLO EW corrections to most of these processes are not yet studied, although they are 
certainly significant in some cases. Therefore, precision theoretical predictions including 
both the NLO QCD and NLO EW corrections for the multiple gauge boson productions 
are necessary. 

The triple gauge boson productions are of particular interest because they are sensitive 
to both the triple and quartic gauge couplings (TGCs and QGCs) and thus related to the 
electroweak symmetry breaking mechanism [3,4], Therefore, the measurements of the 
triple gauge boson productions at hadron colliders can provide rich information on the 
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gauge self-interactions and play an important role in searching for new physics beyond 
the SM. All the triple gauge boson productions at hadron colliders, pp —>• WWZ, ZZZ, 
WWW, WZZ, WW-y, ZZ'y, Z 77 , 777 , W 77 and WZ'y, have been studied in the SM up 
to the QCD NLO [5-13], while only the NLO EW correction to the pp —>• WWZ process has 
been calculated [14]. Therefore, the precision study on the VV'V" (V, V',V" = W or Z) 
productions at hadron colliders with subsequent vector boson decays including the NLO 
QCD + NLO EW corrections is desired, and is listed in the Les Houches 2013 high precision 
wish list [15]. 

The WZZ production at the LHC is sensitive to both the triple WWZ coupling and 
quartic WWZZ coupling and thus relevant for studying anomalous gauge couplings [16,17], 
and this production process with leptonic decays may serve as SM background in searching 
for new physics beyond the SM. The NLO QCD correction to the WZZ production at the 
LHC was calculated in Refs. [7] and [ 8 ], while the NLO EW correction to this process with 
subsequent W- and Z-boson leptonic decays has not been investigated so far. In this paper 
we study the NLO QCD + NLO EW corrections to the WZZ production with subsequent 
vector boson decays at the LHC. The rest of this paper is organized as follows: In section |] 
we provide the details of our calculation strategy. The integrated cross sections and some 
kinematic distributions for the pp —>• WZZ + X process up to the QCD and EW NLO are 
presented and discussed in section ||. Finally, a short summary is given in section [|. 

2. Calculations 

We only take into account the Cabibbo-Kobayashi-Maskawa (CKM) mixing between the 
first two quark generations since the mixing to the third generation is negligible, i.e., the 
CKM matrix is 2 © 1 block-diagonal. The masses of the first two generations of quarks are 
set to zero. In this approximation, the CKM matrix factorizes from all the amplitudes, 
including the tree-level amplitudes for WZZ, WZZ + g, WZZ + 7 , WZZ + q productions 
and the QCD and EW one-loop amplitudes for WZZ production. Therefore, only one 
generic amplitude for each category mentioned above has to be evaluated when convoluting 
the squared matrix elements with the parton distribution functions (PDFs) [18]. We adopt 
the ’t Hooft-Feynman gauge and the four-flavor scheme in the calculations for the LO and 
NLO QCD + NLO EW corrections. 

2.1 LO calculation 

The LO contributions to the W + ZZ and W~ ZZ productions at the LHC come from the 
following partonic processes: 

Qiipi) + Q2(P2) -t W + (p 3 ) + Z(p 4 ) + Z(p 5 ), 

Qi (pi) + <? 2 (P 2 ) ->• w~(p 3 ) + Z(p A ) + Z(p 5 ), (qi=u,c, q 2 = d,s), (2.1) 

respectively. The parton-level cross section for the q 2 qi —> W~ ZZ process in the SM should 
be the same as for the q\q 2 —>• W + ZZ process due to the CP conservation. Therefore, we 
describe the LO and NLO calculations only for the W + ZZ production in this section. 
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The LO Feynman diagrams for the q A q 2 —» W + ZZ partonic process are shown in Fig-0- 
The WWZ TGC is involved in Figs.||(a,b,c,d,e,f) and only Fig.||(g) contains the WWZZ 
QGC. The LO parton-level cross section for q\q 2 —>• W + ZZ is expressed as 
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( 2 . 2 ) 


where the factor ^ arises from the two identical Z-bosons in the final state. The summation 
is taken over the spins of the final state, and the bar over the summation represents 
averaging over the spins and colors of the initial state. -Mq% 0 is the LO Feynman amplitude 
for the qiq 2 —>• W + ZZ partonic process, and dfi 3 is the three-body final state phase space 
element defined as 
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Figure 1: The LO Feynman diagrams for the qiq 2 W + ZZ partonic process, where H and G 
represent the Higgs and charged Goldstone bosons, respectively. 


The final produced W + - and Z-bosons are unstable particles, and we only consider 
their leptonic decay modes in investigating the W + ZZ production. Then the LO cross 
section for the pp —>• W + ZZ —>• l\t 3 Z 3 +X process can be obtained by performing 
the PDF convolution and applying the narrow width approximation as 

°lo(pp -> w + zz + x, s) 
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^ qi \ Fl { x l, /d F )^q 2 \P 2 {x 2 , p F ) + {Pi ^ P 2 ) °Lo{* = x l x 2$) 
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£=e,fj,,T 


£=e,p,T 


(2.4) 


where $ g |p is the quark density in proton, p F is the factorization scale, and \/I and \J~S 
are the center-of-mass system energies of the initial q\q 2 quark pair and colliding protons, 
respectively. 


2.2 NLO QCD corrections 

The NLO QCD correction to the parent process pp —>• W + ZZ + X includes: (1) The QCD 
one-loop virtual corrections to q\q 2 —> W + ZZ] (2) The real gluon emission corrections from 
q i q 2 —> W + ZZ + g: (3) The real light-quark emission corrections from q\g —>• W + ZZ + q 2 
and q 2 g —> W + ZZ + qi, (where q\ = u, c, q 2 = d,s ); (4) The corresponding contributions 
of the PDF counterterms. 

An ultraviolet (UV) and infrared (IR) safe observable requires the exact cancelation 
of the UV and IR singularities to all orders. In this paper, we adopt the dimensional 
regularization scheme in D = 4 — 2e dimensions to isolate both the UV and IR singularities. 
After renormalizing the related quark wave functions ip qi and ip q2 , the QCD one-loop virtual 
correction to the q\q 2 —> W + ZZ partonic process is UV-finite but still contains soft and 
collinear IR singularities. The soft IR singularity is canceled exactly by that in the real 
gluon emission correction from the q±q 2 —>• W + ZZ + g partonic process, while the collinear 
IR singularity is only partially canceled by that in the real gluon emission correction from 
QiQ -2 —> W + ZZ + g and the remaining collinear IR singularity is absorbed by the collinear 
gluon emission parts of the related quark PDF QCD counterterms. The real light-quark 
emission corrections from q\g —> W + ZZ + <72 and q 2 g —>• W + ZZ + qi only contain collinear 
IR singularities, and can be canceled exactly by the collinear quark emission parts of 
the PDF QCD counterterms for q\ and q 2 , respectively. Therefore, A &qq D , Auq^ d and 
A a QCDi defined as 


Aa QCD = J dxidx 2 
+ [ dx\dx2 


z.qiq2,V , ±.qiq2,R 
a QCD “T °QCD 


^ qi \ Pl {xi,p F )^q 2 \p 2 {x2,p F ) + (Pi O P 2 ,X 1 ^ X 2 ) 
5<$>QCB( 9 \ x l,p F )<f>q 2 \P2( x 2, MF) + ^i|Pi {x 1 ,p F )6^p/ 9 \x2,PF) 
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are both UV- and IR-finite separately, where a q(?d , &qcd*J ®qcd an< ^ ®QCD are the 
NLO QCD virtual and real emission corrections from qiq 2 —>• W + ZZ, q\q 2 —> W + ZZ + g, 
qig —>• W + ZZ+q 2 and g 2 g —>• W + ZZ+qi, respectively. The quark PDF QCD counterterm 
5&®p D can be split into two parts: the collinear gluon emission part and collinear light- 

quark emission part, i.e., which can be expressed in the 

MS factorization scheme as 


5$Qp D ' i9 \x,n F ,g R ) 

5 ^\p D ' {q \ x ^ F ^R) 


1 |" C F a s r(l — e) 
e [ 2 t r r(l — 2 e) 
1 ' T F a s T(1 - e) 
e [ 2 vr T (1 - 2 e) 




- 7 - [Pqq(z)]+®q\p(x/z,fl F ), 
dz 

—P qg (z)<f>g\p(x/z,n F ), (2.6) 


where C F = 4/3, T F = 1/2 and is the renormalization scale. The splitting functions 
are given by 


Pqq( z ) ~ 


1 + z 2 
1 - z 


Pqg(z) = z 2 + ( 1 - zf 


(2.7) 


and the [...], prescription is understood as 

z)\ + f{z) = f dzg(z) [f(z) - /(l)]. ( 2 . 8 ) 

Jo 

In our NLO calculation we adopt the two cutoff phase space slicing (TCPSS) technique 
to isolate the IR singularities for the real emission partonic processes [19]. Two cutoffs 
5 S and 5 C are introduced to separate the phase space of the real gluon emission process 
Qi(pi)+q2(p2) -> W + (p 3 )+Z(p A )+Z(p 5 )+g(p 6 ) into soft gluon region (E 6 < 5 s \fi/2), hard 
collinear region (Eq > 8 s V§/2 and min{si6, ^26} < S c s) and hard noncollinear region (Eq > 
8 s V§/2 and min{si6, -§ 26 } > S c s). The real light-quark emission processes, q\ (p\) + <?(p 2 ) —> 
W + (p 3 ) + Z(p 4 ) + Z(p 5 ) + q 2 {pe) and g 2 (pi) + g(p 2 ) ->• W + (p 3 ) + Z(p 4 ) + Z(p 5 ) + qi(p 6 ), 
contain only collinear IR singularities, therefore we only separate their phase space into 
collinear (s 2 6 £ 8 c s) and noncollinear (s 2 6 > 8 c s ) regions. Then we can express the cross 
sections for these real emission partonic processes as 



. - 9152 ,s , ~qiq 2 ,HC , ~qiq 2 ,HC 
°QCD ' ° QCD F °QCD 1 


z.qiq2,R 
a QCD 

~qig,R _ ~qig,C . ~qig,C 
a QCD ~ a QCD "T °QCD > 


±q2g,R _ ~q2g,C , ~q 2 g,C 
a QCD — °QCD "r °QCD > 


(2.9) 


where the superscripts S, HC, HC, C and C stand for soft, hard collinear, hard non¬ 
collinear, collinear and noncollinear, respectively. The cross sections over the soft and 
(hard) collinear regions contain only the soft and collinear IR singularities separately, while 
the cross sections over the (hard) noncollinear regions are IR-finite. We have checked nu¬ 
merically the cutoff independence of the total cross sections for these subprocesses by 
setting 5 C = 5 S /50 and varying 5 S from ICR 5 to 10~ 3 . 

Finally, we get the full NLO QCD correction to the parent process pp —>• W + ZZ + X 
as 


Ao-qcd 


q 2 =d,s 

£ [ 

qi=u,c 


A * 7 QCD + A CTqcd + A Vqcd 


( 2 . 10 ) 
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We employ the modified FeynArts 3.7 package [20] to generate Feynman diagrams and 
their corresponding amplitudes. The reduction of output amplitudes are implemented by 
the FormCalc 7.3 package [21]. In our numerical calculation, the tensor 4-point integrals 
with rank n > 3 may induce a serious unstable problem. We adopt our developed codes 
based on the LoopTools 2.8 package [22] to calculate the scaler and tensor integrals, which 
can switch to the quadruple precision arithmetic automatically if necessary. 

In Ref. [7] the authors took the same input parameters as in Ref. [8] and made a 
numerical comparison for the pp —>• W + ZZ + X process at the \/S = 14 TeV LHC. We 
follow their input parameters and settings, and present in Table [l] a comparison between 
our numerical results for the LO and NLO QCD corrected integrated cross sections and 
the corresponding ones provided in Ref. [7] and Ref. [8]. It shows that all these numerical 
results are in good agreement with each other within the Monte Carlo errors. 


pp —*• W + ZZ + X 

ctlo [fb] 

&NLO [fb] 

ours 

20.3(1) 

39.7(2) 

Ref. [7] 

20.24(3) 

39.86(7) 

Ref. [8] 

20.0(1) 

39.7(2) 


Table 1: Comparison of our numerical results for the LO and NLO QCD corrected integrated 
cross sections with the corresponding ones in previous works [7,8]. All input parameters and settings 
are taken from Ref. [7] with pp = pp = 3Mz- 


2.3 NLO EW corrections 


The NLO EW correction to the parent process pp —» W + ZZ+X includes: (1) The EW one- 
loop virtual corrections to q\q 2 —>• W + ZZ] (2) The real photon emission corrections from 
q\ q 2 —>• W + ZZ + 7 ; (3) The contributions of the photon-induced partonic processes 517 —>• 
W + ZZ + (/2 and 527 —>• W + ZZ + qi, (where q\ = u,c, <72 = d,s)\ (4) The corresponding 
contributions of the PDF counterterms. 

The UV divergences in the EW one-loop virtual corrections can be removed by the 
renormalization procedure. We take the definitions for the relevant EW renormalization 
constants same as in Ref. [23]. We adopt the on-nrass-shell scheme to renormalize the 
masses and wave functions of related particles as used in the QCD correction. The expres¬ 
sions for the relevant renormalization constants and the unrenormalized EW self-energies 
can be found in Ref. [23]. In our calculation, we adopt a mixed scheme to deal with the 
EW couplings. All the EW couplings in the tree-level diagrams for < 71(72 —>• W + ZZ are 
fixed in the G^-scheme, i.e., a = ac^, while the extra EW couplings appeared in the EW 
one-loop diagrams for < 71(72 W + ZZ and in the real photon emission and photon-induced 
subprocesses are fixed in the a( 0 )-scheme, i.e., a = a( 0 ) 1 . a( 0 ) is the fine structure 
constant in the Thomson limit and is given by 


«G„ = —-GpMw sin 2 9w 
7r 


( 2 . 11 ) 


1 The fine structure constant in the PDF EW counterterms (see Eq.(2.15)) is also fixed as a = a(0). 
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where sin 2 9\y = 1 — Af /M J-. Then the LO cross section and NLO EW correction are of 
O(ctQ^) and 0(o.q a( 0 )), respectively, and correspondingly the electric charge renormal¬ 
ization constant should be given in the Gp-scherne as 

sz^ = §zf°'> - I A r = —SZaa - \ tan 9 w 5Z za - i Ar, ( 2 . 12 ) 


where SZe^ is the electric charge renormalization constant in the a(0)-scheme, and Ar 
[23, 24] corresponds to the subtraction of the logarithmic divergence contributed by the 
light quarks to SZ AA which was absorbed by olq . 

The quark PDF EW counterterm <54>^jp also contains two parts: the collinear photon 

emission part and collinear light-quark emission part, i.e., . 

In the DIS factorization scheme used for NLO EW corrections, these two collinear parts 
are expressed as [14,25] 


S^^\ x ,h f ,hr) 

S^\^ q \ Xj HF,n R ) 



1 r(l - e) 

e T (1 - 2 e) 

f 1 r(l-e) 

l 6 T (1 — 2 e) 



4 ^/?. 


[P qq (z)] + -C° IS (z) 
Pq'r(z) - C° IS (z) | 


( 2 - 


where Q q is the electric charge carried by the initial quark q, and P qi (z) = P qg {z). The 
DIS factorization scheme is specified by [14,25] 


fl*) = 


Pqq(z) In 


1 - 2 


+ 


9 + 5z 


C?J s (z) = P qi (z) In 


1 - z 


•<?7 


- 8z 2 + 8z - 1. 


(2.14) 


Analogous to the QCD correction described in Sec. |2.2| , we get the full NLO EW correction 
to the parent process pp —>• W + ZZ + X as 

q2=d,s 


Actew = E [ ^ a< EW 

qi=u,c 


(2.15) 


where the analytic expressions for Adp 1 ^, Acrpjy and Aa q ^ v can be obtained from Eqs.(2.5) 
and ( |2.9D by doing the replacements of QCD —> EW and g — 7 . Each of Acr^ff, A 
and Aap(y ( q\ = u, c, q 2 = d, s ) is UV- and IR-hnite, therefore the full NLO EW correction 
A ctew is an UV- and IR-safe variable. 


2.4 CKM matrix dependence 

The W + ZZ production at the LHC up to the QCD and EW NLO involves the following 
topologies: 


0 —> W + ZZ + q\q 2 
0 —> W + ZZ + qiq 2 + g 
0 —> W + ZZ + qiq 2 + 7 


(tree, QCD and EW loop) 
(tree), 

(tree), 


(2.16) 
















where q\ = u,c and q 2 = d,s. Each of these topologies contains only one charged current 
quark chain , 2 therefore the CKM matrix element can factorize from the amplitudes as 


■ M . tree, loop (0 — 5 - W + ZZ + — V qi q 2 X A 7 tree, loop (0 — > W + ZZ + ud ), 

-Adtree(0 —> W + ZZ + q±q2 + g) = V qi q 2 X -^tree(0 —* W + ZZ + ud + g), 

M tree ( 0 -> 1E+ZZ + q x q 2 + 7) = Vq iq2 X Mree(0 ->• PE+ZZ + ud + 7 ), (2.17) 

since m u = m c = = m s = 0, where M and M are the amplitudes obtained with and 

without CKM matrix. 

For convenience we introduce the following CKM matrix dependent structure functions: 


q 2 =d,s 

Fqq{xi,X2,H F ) = ^ \Vqiq 2 \ 2 ®qi |Pi (®1, /d F )^q 2 \P 2 {x 2 , Vf) + (-Pi ^ -?2, X X -B- X 2 ) 
qi=u,c 
q2=d,s 

J r qg(x 1 ,X 2 ,IJ,F) = ®qi\Pi ( x l, Mf)$ 9 |P 2 (x 2 , Vf) + (A ^ ^2, X\ -B X 2 ) 

qi=u,c 
q 2 =d,s 

Fqgixi, X 2 , Hf) = ^ \ V qiq 2 \ 2 $ 5 |Pi ( x l, ^ F )^q 2 \P 2 ( X 2, Mf) + (A ^2, Xl B- X 2 ) 

qi=u,c 


d~q^ (lC 1 , %2 7 /^f) 
J r q 1 {xi,X2, Hf) 


J 7 qg(xi,X2, fiF ) 

Fqg(x 1, X2, [Ip] 


ff->7 


(2.18) 


q—t'y 

Then the full NLO QCD and EW corrections to the parent process pp —> 1E + ZZ + X can 
be simply rewritten as 


Act qcd = J dxidx 2 
Act ew = / dx\dx2 


Fan a, 


■u d,V 
qq 1 ^qcd 


+ cr, 


ud,R \ 
QCD J 


-r- ~ug,R _■ 

+ -F qg®QCD T •'aa&t 


dg,R 
qg u qcd 


i c -pQCD 

+ 0-F q g <r LO 




wd 


/ *ud,V zud,R\ EW _ ™ 

■'qq I ^ I + -F q-y cr EW + J ~q'Y Cr EW + °-F qq a LO 


(2.19) 


where 5F^ D and SF^ g w are the QCD and EW counterterms of respectively, expressed 


99 
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&7v 


QCD 


q 2 =d,s 
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E ib 


9192 I 


( 1 $ 


QCD 


j^(®l» /^f)^ 52 |p 2 (x 2 ,hf) + $ qi I Pi (® 1 , (® 2 , Mf) 


91 


qi=u,c 


+ (Pi B P 2 ,®1 -H- X 2 ) 


r jrEW r -pQCD 

0J ~qq 0J ~qq 


QCD^EW 


( 2 . 20 ) 


^ ud zud,V zud,R zug,R zdg.R, zud,V zud,R ^ u'y^R ^dj,R 
and cr L Q, ctq CD , (Tqcd > a QCDi a QCDi a EW > ^FVF > a EW ’ ^FVF are tae corresponding 

partonic cross sections by setting Vckm = 7. From Eq.( 2 . 1 £) we see clearly that the 


2 We do not consider the closed quark loop, because the CKM matrix in it drops out after the summation 

over quark flavors. 



















CKM matrix factorizes from the amplitudes and is absorbed by the structure functions. 
Therefore, only the amplitudes for one generation of quarks have to be evaluated in PDF 
convolution 3 . 


3. Numerical results and discussion 


3.1 Input parameters 

The SM input parameters are taken as [26]: 


M w = 80.385 GeV, M z = 91.1876 GeV, M H = 126 GeV, m t = 173.5 GeV, 

a(0) = 1/137.036, a s (m z ) = 0.119. (3.1) 


G^ = 1.16638 x 10" 5 GeV" 2 



f v ud 

V us 

v ub \ 

/ 0.97425 

0.22547 

0 ^ 

Vckm = 

v cd 

v cs 

V cb 

= -0.22547 

0.97425 

0 


\v td 

V ts 

Vtb 

l 0 

0 

V 


The masses of all leptons and quarks but the top quark are set to zero. The CKM matrix 
elements are taken as 


(3.2) 


We set the factorization and renormalization scales being equal and choose the central scale 
as ho = + M Z - The two cutoffs are set as 5 S = 10 —3 , 5 C = 2 x 10^ 5 in NLO QCD 

calculation, and 5 S = 10 , 6 C = 2 x 10~ 6 in NLO EW calculation, respectively. 

In the calculation of the NLO QCD corrections, we adopt the NLO NNPDF2.3QED 
PDFs [27] with the MS factorization scheme. The strong coupling constant a s is renor¬ 
malized in the MS scheme with five active flavors, and its running value a s (/u) is taken 
from the PDF set. While for the NLO EW corrections we use the NLO NNPDF2.3QED 
PDFs with the DIS factorization scheme [25]. 


3.2 Integrated cross sections 

The NLO QCD and EW relative corrections are defined as 

Aoqcd + (oo - Tlo) 


Sqcd = 


&LO 


5ew = 


Act EW 
Oo 


(3.3) 


where Actqcd and Aoew (see Eqs.( |2.ld ) and ( 2.15 )) are evaluated with NLO PDFs, and 
c teo and &o are LO cross sections calculated with LO and NLO PDFs, respectively. The 
numerator Act qcd + (co — &lo) represents the full NLO QCD correction that includes 
all the NLO QCD contributions from both the dynamic matrix element and PDFs. To 
cancel the QCD contribution from NLO PDFs to the NLO EW correction ActeWi we 
normalize the NLO EW relative correction to oo- hr this normalization, the NLO EW 
relative correction 5ew is practically independent of the PDF set. 


3 In the case of Vckm = I and m u = m c = mj = m s = 0, the amplitudes for the first quark generation 
are the same as the corresponding ones for the second quark generation. In this paper we compute the 


related amplitudes and partonic cross sections only for the first generation of quarks (see Eqs.(2.16) and 

©)■ 
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Ideally, the NLO QCD+EW correction should be calculated by applying the multi-jet 
merging approach, which would imply the calculation of NLO QCD and EW corrections 
to W±ZZ + jet final states and thus is beyond the present scope. Therefore, in this paper 
we calculate the combined NLO QCD+EW correction by using the naive product [28] 

&nlo = Tlo (1 + Snlo) 

= olo (1 + &qcd) (1 + Sew) , (3.4) 

where (Jnlo is the NLO QCD+EW corrected cross section and 5nlo is the NLO QCD+EW 
relative correction. 

In order to keep the convergence of the perturbative QCD description of the pp —>• 
W^ZZ + X processes, we may impose a tight jet veto which can heavily suppress the large 
QCD correction. We call the event selection scheme with jet veto condition of pr,jet < 
P^Tjet = 50 GeV as the exclusive scheme (scheme-II) and that without any jet veto as the 
inclusive scheme (scheme-I). In Tableland Table ||, we present the LO and NLO QCD+EW 
corrected cross sections for pp —>• W±ZZ + X in both the inclusive and exclusive event 
selection schemes at the 13 and 14 TeV LHC, respectively 4 . The EW corrections from 
the quark-antiquark and photon-induced channels are shown separately, since they can (in 
principle) be distinguished by their final states. 

From Table [2] and Table |3] we can see that the photon-induced channels have surpris¬ 
ingly large impact on the NLO EW correction. The NLO EW correction can be heavily 
suppressed by applying a jet veto. For example, the photon-induced EW relative correc¬ 
tion 5 q ^j w is 16.23% in the inclusive event collection scheme, but is reduced to 1.73% after 
applying pr,jet < 50 GeV on the final jet, for the pp —» W + ZZ + X process at the 13 TeV 
LHC. It implies that the theoretical uncertainty from the photon PDF can be reduced by 
adopting the jet-veto event selection scheme. In our calculation we find that the real light- 
quark emission correction is the largest NLO QCD contribution and amounts to 56.3% 
of the full NLO QCD correction for the W + ZZ production at the 14 TeV LHC. In both 
two event collection schemes, the quark-antiquark and photon-induced EW corrections are 
negative and positive, respectively. The full NLO EW relative correction to the W + ZZ 
production in the inclusive event collection scheme at the 14 TeV LHC can reach about 
9.67%. 

Now we turn to the fractorization/renormalization scale dependence of the LO and 
NLO corrected integrated cross sections. The factorization scale pf affects both the LO 
and NLO corrected cross sections via the factorization procedure due to the pf dependence 
of the parton densities. The renormalization scale pp occurs in higher order perturbative 
calculation via the renormalization procedure and strongly affects the QCD correction. For 
simplicity, we set pf = pp = p in our calculation. In Figs.||(a) and (b), we depict the 
LO, NLO QCD and QCD+EW corrected cross sections (o+o> &QCD and ctnlo ) f° r the 
W + ZZ production at the 14 TeV LHC as functions of p by adopting the inclusive and 
exclusive event selection schemes, respectively. The corresponding NLO QCD+EW relative 

4 In Table || and Table 1 Ao-qcd should be understood as the full NLO QCD correction, i.e., A uqcd = 
ctloSqcd- 
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process 

pp W+ZZ + X 

pp W 

ZZ + X 

scheme-I 

scheme-II 

scheme-I 

scheme-II 

0+0 [fb] 

17.078(5) 

17.078(5) 

8.625(2) 

8.625(2) 

CTNLO [fb] 

42.71(5) 

20.68(3) 

24.19(2) 

11.39(2) 

AtTQCD [fb] 

22.06(4) 

4.79(3) 

12.82(2) 

3.29(2) 

Aa< EW if b ] 

-1.412(6) 

-1.412(6) 

-0.710(3) 

-0.710(3) 

Acr EW if b \ 

3.222(2) 

0.3443(4) 

2.061(1) 

0.2335(3) 

$QCD [%] 

129.2 

28.0 

148.6 

38.1 

b q E q w [%] 

-7.11 

-7.11 

-6.74 

-6.74 

b< EW [%] 

16.23 

1.73 

19.57 

2.22 

5nlo [%] 

150.1 

21.1 

180.5 

32.0 


Table 2: The LO, NLO QCD+EW corrected integrated cross sections for pp —> W+ZZ + X and 
the corresponding NLO corrections at the = 13 TeV LHC. 


process 

pp W+ZZ + A' 

pp —>■ W 

ZZ + X 

scheme-I 

scheme-II 

scheme-I 

scheme-II 

O+O [fb] 

19.126(5) 

19.126(5) 

9.816(3) 

9.816(3) 

CTNLO [fb] 

48.96(4) 

23.05(3) 

28.19(2) 

12.87(2) 

a cjqcd [fb] 

25.52(4) 

5.25(3) 

15.05(2) 

3.67(2) 

Aa E:W [fb] 

-1.601(6) 

-1.601(6) 

-0.819(3) 

-0.819(3) 

A(jq EW [fb] 

3.753(2) 

0.3952(5) 

2.424(1) 

0.2702(3) 

$QCD [%] 

133.4 

27.4 

153.3 

37.4 

S q E q w [%] 

-7.20 

-7.20 

-6.83 

-6.83 

bEW [%] 

16.87 

1.78 

20.22 

2.25 

Snlo [%] 

156.0 

20.5 

187.2 

31.1 


Table 3: The LO, NLO QCD+EW corrected integrated cross sections for pp —> W+ZZ + X and 
the corresponding NLO corrections at the \fS = 14 TeV LHC. 


corrections are shown in the lower plots. To estimate the theoretical uncertainty from the 
factorization/renormalization scale quantitatively, we introduce the scale uncertainty as 


max {cr(/i)|0.2/xo < /U < 5^o} — min{<r(/it)|0.2//o < ^ < 5/+)} 

^(Mo) 


(3.5) 


From Figsj| (a,b) we obtain p = 1.6% at the LO, and r/ = 10.8% and 7.5% at the QCD 
NLO in the inclusive and exclusive event collection schemes respectively. We can see that 
the scale uncertainty at the LO is much less than at the QCD NLO, because the strong 
coupling a s is not involved in the LO matrix elements. When we include both the NLO 
QCD and NLO EW corrections into consideration, the scale uncertainties are 7.9% and 
8.4% in the inclusive and exclusive event selection schemes, respectively. We can conclude 
that the scale uncertainty at the QCD+EW NLO mainly comes from the QCD correction. 


From the numerical results given in this subsection, we can draw the following four 
conclusions: (1) The NLO QCD correction is very large and ruins the convergence of the 
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Figure 2: The fractorization/renormalization scale dependence of <tlo (solid), ctqcd (dotted) and 
&nlo (dashed) for the W + ZZ production at the \JS = 14 TeV LHC in the (a) inclusive and (b) 
exclusive event selection schemes. 


perturbative QCD description in the inclusive event selection scheme. (2) The NLO QCD 
correction mainly comes from the real jet radiation subprocesses, and therefore can be 
heavily reduced by adopting the exclusive event selection scheme with a jet veto. But a 
tight jet veto would introduce a new source of theoretical uncertainties from various other 
processes. (3) Due to the small scale dependence of PDFs in the Feynman-x region, the 
scale dependence of the LO cross section is not apparent. But the LO scale uncertainty does 
not give a good estimate of the higher order QCD contribution. (4) Compared to the NLO 
QCD correction, the NLO EW correction is insensitive to the factorization/renormalization 
scale. 

3.3 Kinematic distributions 

In this subsection we present some kinematic distributions of final produced particles for 
the W + ZZ production at the 14 TeV LHC. In order to take into account the off-shell con¬ 
tributions and spin correlations from the W + - and Z-boson leptonic decays, we transform 
the differential cross sections into Les Houches event files [29] and use MadSpin [30] to 
obtain events after the vector boson decays. For each kinematic variable x considered in 
the following, we provide the LO, NLO QCD and QCD+EW corrected distributions, i.e., 
d&Lo/dx, doQCD/dx and da^Lo/dx. 

We present the invariant mass distributions of W + ZZ system by adopting the inclusive 
and exclusive event collection schemes in Figs.^(a) and (b), respectively. The corresponding 
relative corrections are depicted in the nether plots. The figures show that the W + ZZ 
invariant mass distributions in both event selection schemes reach their maxima at the 
position of M w + zz ~ 440 GeV, and the inclusive NLO QCD correction enhances the LO 
Mw+zz distribution significantly. From Fig||(a) we see that the NLO EW correction in the 
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inclusive event selection scheme is positive but very small compared with the corresponding 
NLO QCD correction in the region of M w + zz > 400 GeV. From Fig.||(b) we find that 
in the exclusive event selection scheme the NLO QCD correction is heavily reduced and 
the NLO EW correction becomes negative because of the jet veto in the exclusive event 
selection scheme. In the plotted M w +zz region, the NLO QCD+EW relative corrections in 
the inclusive and exclusive event selection schemes range from 123% to 170% and from —1% 
to 38%, respectively. The NLO QCD+EW relative correction to the M w + zz distribution 
in the inclusive event selection scheme is more dependent on the phase space than in the 
exclusive event selection scheme, because the NLO QCD correction from the gluon-induced 
channels is the dominant contribution at the NLO and is more closely related to the phase 
space. We conclude that the real jet radiation would induce large NLO contribution in the 
inclusive event collection scheme, and we can keep the convergence of the perturbative QCD 
description and obtain moderate NLO QCD+EW correction by adopting the exclusive 
event selection scheme. However, a jet veto as tight as pr,jet < 50 GeV introduces a 
new source of theoretical uncertainties from various processes. Such uncertainties could be 
analyzed by dedicated jet veto resunnnations. 




Figure 3: The W + ZZ invariant mass distributions daLo/dM w+zz (solid), daQc d/ dM w + zz 
(dotted), d<TNLo/dM w + zz (dashed) and the corresponding relative corrections for pp W + ZZ + 
X at the y/S = 14 TeV LHC in the (a) inclusive and (b) exclusive event selection schemes. 


The rapidity distributions of W + ZZ system by adopting the inclusive and exclusive 
event collection schemes are shown in Figs.[|(a) and (b) separately. The corresponding 
relative corrections are plotted in the lower panels. From the figures we can see that the 
NLO QCD relative corrections in the inclusive and exclusive event collection schemes are 
about 165% and 33%, respectively, at the position of Vw+zz = 0- The NLO EW correction 
enhances the LO W + ZZ rapidity distribution a little bit in the inclusive event selection 
scheme, but reduces the LO Vw+zz distribution in the exclusive event selection scheme. In 
the inclusive event selection scheme the NLO QCD correction is much larger than the NLO 
EW correction. However, the NLO QCD correction in the exclusive event collection scheme 
is heavily reduced, particularly in the region of \yw+zz\ < 1- As we know, the events with 
large W + ZZ transverse momentum, i.e., Pt,w+zz > 50 GeV, tend to be produced more 
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centrally, i.e., Uw+zz 0 , and will be excluded in the exclusive event selection scheme due 
to the jet veto and the conservation of transverse momentum. That’s why the NLO QCD 
correction in the exclusive event selection scheme is small in the range of \yw+zz\ < 1- 
The NLO QCD+EW correction in the exclusive event selection scheme suppresses the LO 
Uw+zz distribution in the vicinity of yw+zz = 0- 
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Figure 4: The W + ZZ rapidity distributions daLo/dyw+zz (solid), daQCD/dy w + Z z (dotted), 
daNLo/dyw+zz (dashed) and the corresponding relative corrections for pp —> W + ZZ + X at the 
y/S = 14 TeY LHC in the (a) inclusive and (b) exclusive event selection schemes. 


The Z-pair invariant mass distributions by adopting the inclusive and exclusive event 
collection criteria are plotted in Figs.|5](a) and (b), respectively. Form the figures we see 
that the NLO QCD and EW corrections do not distort the line shape of the LO Mzz 
distribution, and both the LO and NLO corrected Mzz distributions reach their maxima 
at Mzz ~ 210 GeV. The NLO EW correction slightly enhances and suppresses the LO 
Mzz distribution in the inclusive and exclusive event collection schemes, respectively, in 
the plotted Mzz region. The NLO QCD relative correction can exceed 138% when Mzz > 
300 GeV in the inclusive event collection scheme, and is less than 36% in the whole plotted 
Mzz region in the exclusive event collection scheme. It indicates that the NLO QCD 
correction to the Mzz distribution in the inclusive event collection scheme could be very 
large and destroy the perturbative description, and the perturbative convergence can be 
improved by applying a tight jet veto in the exclusive event collection scheme. 

Concerning the sequential leptonic decays of the final W + and Z-bosons, we study 
the pp —» W + ZZ —>• £3 + X process by adopting the improved narrow width 

approximation. The branch ratios for the W- and Z-boson leptonic decay modes are 
obtained by using the MadSpin program. We define the final lepton with the largest 
transverse momentum among all leptons as the leading lepton and that with the second 
largest transverse momentum as the next-to-leading lepton. In Figs|6|(a,b), Figsj^(a,b) 
and Figs.||(a,b) we display the distributions of the leading lepton, next-to-leading lepton 
and missing transverse momenta in the inclusive and exclusive event collection schemes 
separately. The corresponding relative corrections are drawn in the nether panels. We see 
from these figures that all the transverse momentum distributions of the leading lepton, 
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Figure 5: The Z -pair invariant mass distributions daLo/dMzz (solid), dagco/dMzz (dotted), 
d(JNLo/dMzz (dashed) and the corresponding relative corrections for pp —> W + ZZ + X at the 
y/S = 14 TeV LHC in the (a) inclusive and (b) exclusive event selection schemes. 


next-to-leading lepton and missing energy have similar behavior. Both the NLO QCD and 
NLO EW corrections enhance the LO pp distributions in the whole plotted pp region in the 
inclusive event selection scheme, but suppress the LO pp distributions in high pp region 
in the exclusive event selection scheme. In the inclusive event selection scheme, the NLO 
QCD relative correction increases with the increment of pp in the region of pp > 50 GeV. 
For example, we can read out from Fig|](a) that the NLO QCD relative correction in the 
inclusive event selection scheme increases from 82% to 235% as the increment of p!f lcp from 
50 GeV to 250 GeV. The large NLO QCD correction in high pp region in the inclusive 
event selection scheme is dominated by the gluon-induced channels [7,14]. In the exclusive 
event selection scheme, the NLO EW correction is very small and negligible compared with 
the corresponding QCD correction in the whole plotted pp region. The NLO QCD relative 
correction increases in low pp region and then drops down in high pp region since most of 
the high-pT events are rejected by the jet veto. 

In above calculation we neglect the NLO EW correction to the pp —>• W + ZZ —» 
e+v tl i+qe+i j + x process from the real photon radiation off the final-state charged 
leptons. As we know, photon radiation off a final-state charged lepton reduces the lepton 
momentum and thus shifts many events out of the acceptance window as well as modifies 
the kinematic distributions of final products. This correction can be significant due to the 
mass-singular logarithms log — induced by the small lepton mass. We will include this 
contribution in leading logarithmic accuracy in future improved calculation. 

4. Summary 

In this paper, we investigate the NLO QCD + NLO EW corrections to the W±ZZ pro¬ 
ductions with subsequent W±- and Z-boson leptonic decays at the 14 TeV LHC, by adopt¬ 
ing the MadSpin method which preserves both spin correlation and finite width effects 
to a very good accuracy. The NLO QCD+EW corrected integrated cross section and 
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Figure 6: The leading lepton transverse momentum distributions do lo t dp^Z lep (solid), 
d<JQCD/dp lep (dotted), doNLo/dp ^ lep (dashed) and the corresponding relative corrections for 
pp —> W + ZZ —> at the \[S = 14 TeV LHC in the (a) inclusive and (b) exclusive 

event selection schemes. 




Figure 7: The next-to-leading lepton transverse momentum distributions doLo/dpi^ L lep (solid), 
dcrQCD/dp^ L lep (dotted), do^Lo/dp^ 1 " lep (dashed) and the corresponding relative corrections 
for pp —> W + ZZ —> + X at the \/~S = 14 TeV LHC in the (a) inclusive and (b) 

exclusive event selection schemes. 


some kinematic distributions are studied. Our results demonstrate that the NLO QCD 
and NLO EW corrections are all significant, and modify the LO integrated cross section 
and some kinematic distributions obviously. In the jet-veto event selection scheme with 
P^Tjet = 50 GeV, the NLO QCD+EW relative corrections to the integrated cross section 
are 20.5% and 31.1%, while the genuine NLO EW relative corrections are —5.42% and 
—4.58%, for the W + ZZ and W~ZZ productions, respectively. We also investigate the 
factorization/renormalization scale dependence of the integrated cross section, and find 
that the scale uncertainty is underestimated at the LO due to the absence of the strong 
coupling a s in the LO matrix elements. 
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Figure 8: The missing transverse momentum distributions doLo/dpx lss (solid), daQCD / dp™ lss 
(dotted), daNLo/dp™ lss (dashed) and the corresponding relative corrections for pp — > W + ZZ —► 
+ X at the y/S = 14 TeV LHC in the (a) inclusive and (b) exclusive event selection 

schemes. 
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